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Abstract 

We have searched for mixmg m the D^-D^ system by measuring the decay-time distribution 
of — > K~^n~ decays. The analysis uses 90 fb~^ of data collected by the Belle detector at the 
KEKB e+e" collider. We fit the decay-time distribution for the mixing parameters x' and y' and 
also for the parameter Rj-,, which is the ratio of the rate for the doubly-Cabibbo-suppressed decay 
^ K'^tt" to that for the Cabibbo-favored decay D^^K^tt'^. We do these fits both assuming 
CP conservation and allowing for CP violation. We use a frequentist method to obtain a 95% 
C.L. region in the x'^-y' plane. Assuming no mixing, we measure i?^, = (0.381 ± 0.017 ^o g^e)^- 

PACS numbers: 12.15.Ff, 13.25.Ft, ll.30.Er 
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The phenomenon of mixing among quark flavors has been observed in the K^-K^ and 
B^-B^ systems but not yet in the D^-D^ system. The rate for D^-D^ mixing within the 
Standard Model (SM) is small, typically well below experimental upper limits yj. Obser- 
vation of mixing much larger than this expectation could indicate new physics, such as a 
AC = 2 interaction. Such nonstandard processes may also give rise to CP- violating effects. 

In this paper we present a search for D^-D^ mixing and CP violation {CPV) in mixing 
with greater sensitivity than that of previous searches. The data sample consists of 90 fb~^ 
recorded by the Belle experiment at the KEKB asymmetric e~^e~ collider The Belle 
detector |2| consists of a three-layer silicon vertex detector (SVD), a 50-layer central drift 
chamber (CDC), an array of aerogel threshold Cherenkov counters (ACC), time-of-fiight 
scintillation counters (TOF), and an electromagnetic calorimeter comprised of CsI(Tl) crys- 
tals. These detectors are located within a solenoid coil providing a 1.5 T magnetic field. An 
iron fiux-return outside the coil is instrumented to identify muons and mesons. 

The dominant two-body decay of the is the Cabibbo-favored (CF) decay — >■ 
K~7i^ 01 • We search for mixing by reconstructing the "wrong-sign" decay — > K~^7i~, 
which would arise from a mixing to and subsequently decaying via ^ K^tt~ . 
The flavor of the D is identified by requiring that it originate from D*^ D^tc^ or 



D* 



D n and noting the charge of the accompanying pion. In addition to arising 



via mixing, K^7[~ can also occur via a doubly-Cabibbo-suppressed (DCS) amplitude. 
The two processes can be distinguished via the decay-time distribution. This method has 
been used by FNAL E791 [Si], CLEO and BaBar i2] to search for mixing and measure 
or constrain the DCS decay rate. 

The parameters used to characterize mixing are x = Am/T and y = Ar/(2r), where 
Am and AF are the differences in mass and decay width between the two D^-D^ mass 



eigenstates, and F is the mean decay width. For \y\ 
time distribution for —>■ K~^7i~ can be expressed as 



^ 1 and negligible CPV, the decay 
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where Rjy is the ratio of DCS to CF decay rates, x' = x cos 5 + ysmS, y' = y cos 5 — x sin 5, 
and 5 is the strong phase difference between the DCS and CF amplitudes. The first term in 
brackets is due to the DCS amplitude, the last term is due to mixing, and the middle term 
is due to interference between the two processes. The time-integrated rate for — >■ iiT+Tr" 
relative to that for D'^ ^ K~tt~^ is Rr 



[x'^ + y'^)/2 



To allow for CP violation, we follow Ref. 0| and apply Eq. (0) to and decays 



separately. This results in six observables: {-R/), x'"*"^, y'^} for and {Rq, 



X 

?+ 



} for 
and 



. CP violation is parametrized by the asymmetries Aj-^ = (i?J — i?^)/(i?J -|- Rj^ 

(x'^^ + y'^"^) /2. The asymmetry char- 
characterizes CPV in D^-D^ mixing. 



Am = (Rli - Rm)/{Rm + R~m). where R 
acterizes CPV in the DCS decay amplitude, and A^ 
The observables are related to x' and y' via 



X 




X cos ± sin ( 
{y' cos (j)^ x' sin < 
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where is a weak phase and characterizes CPV occurring via interference between mixed 
and unmixed decay amphtudes. Note that x'^, y'^ are unchanged by the transformation 
a;'— x', y'^—y', and 0— s>0 + 7r; thus for definiteness we restrict to the range |0| < n/2. 

We select — > K^n^ decays by requiring two oppositely-charged tracks, with at least 
four SVD hits, that satisfy K and vr identification criteria. These criteria are £^>0.6 and 

> 0.4, where C is the relative likelihood for a track to be a or tt based on dE/dx 
information in the CDC and the responses of the TOF and ACC systems. These criteria 
have efficiencies of 88.0% and 88.5% and tt/K misidentification rates of 8.5% and 8.8%, 
respectively. We combine the candidate with a low-momentum pion (TTgi^^) to form a 
D* + ^D°7r+ candidate. Candidates in which the charge of vr^i^^ is opposite (equal to) that 
of the are referred to as "right-sign" or RS ("wrong-sign" or WS) decays. To reject WS 
background from —>■ K~tt'^ in which the K is misidentified as vr and the vr is misidentified 
as K, we recalculate ruj^^ with the K and vr assignments swapped and reject events with 
I ^K-n~ ^Do\ < 28 MeV/c^. To eliminate D*'s from B decays, we require p^,, > 2.5 GeV/c, 
where p^,* is evaluated in the e+e~ center-of-mass frame. 

The vertex is obtained by fitting the daughter K/n tracks. The D* vertex is taken as 
the intersection of the momentum vector with the interaction profile region. We require 
a good for each vertex fit. The momentum of tt^i^^ is refitted with the constraint that it 
originates from the D* vertex. The D° decay time is calculated as (^^jo/p^jo) x m^o, where 
i^o is the distance between the and D* vertices projected onto the pj^o direction. The 
decay-time resolution is typically 0.2 ps. 

We measure Rj^, x'^, and y' of Eq. ([Q) via an unbinned maximum likelihood fit to the WS 
decay-time distribution. The likelihood function consists of probability density functions 
(pdf's) for signal and several backgrounds. The pdf's depend on the decay time, the mass 
rrij^^, and the kinetic energy released Q = m^^^^ — m^^— m^. The latter equals only 

slow 

5.85 MeV/c^ for D*^ ^ D^n^^^^ Ktitc^^^ decays, which is near the threshold. 

The signal pdf for event i is smeared by a resolution function R- = (1 — f^^^^) G(t^ — 
t', (Tj j ; fi, S) + /^g^ji G{t^ — t', (Tj j ; /i, S'^.j^ii), where the G's are Gaussians with common mean 
/X and standard deviations (5* x J and {S^.^^^ ^ "^t,*)' ^.nd ^ is the uncertainty in decay 
time t for event i. The parameters /ta^jj, /x, and scaling factors S and S'^a^ji are determined 
from data. The background pdf's are smeared by similar resolution functions (see below). 
To check the resolution function, we fit the RS sample in the same manner as the WS sample 
except that the signal pdf has a purely exponential time dependence. 

There are four backgrounds to the WS sample: (a) random vr background, in which a 
random ti^ is paired with a K^n^ decay (the pdf is peaked in rrij^^ but broad in Q); 

(h) D*^ D^TT^ followed by decaying to > 3-body final states (the pdf is broad in rrij^^ 
and broad but enhanced in Q); (c) / decays; and (d) combinatorial. We determine 
the level of each background by performing a two-dimensional fit to the ttlj^^- Q distribution. 
When fitting the RS sample, the ttlj^^ and Q means and widths for signal are fioated; when 
fitting the WS sample, these means and widths are fixed to the values obtained from the RS 
fit. Also for the WS fit, the ratio of D^/Df to > 3-body backgrounds is fixed to the value 
obtained from Monte Carlo (MC) simulation. The RS fit finds 227 721 ± 497 D°^K-7i+ 
decays, and the WS fit finds 845 ± 40 D° ^ K+n- decays. The ratio i?ws = r(D° ^ 
K^7t-)/T{D^-^K-n^) = (0.371 ±0.018)% (statistical errors only). The ratio of WS signal 
to background is 0.9; the latter is mostly random tt (59%) and combinatorial (36%). The 
WS rrij^^ and Q distributions are shown in Fig. ^ along with the fit projections. 

To fit the decay-time distributions of RS and WS samples, we consider the 4a region 
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FIG. 1: Distributions of (a) WS ruj^^ with | Q - 5.9 MeV | < 0.6 MeV; and (b) WS Q with 
' ^Do \ <20 MeV/c^. Superimposed on the data (points with error bars) are projections of 
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Q fit. 



\mj^^ — m^o| < 22 MeV/c^ and \ Q — 5.9 MeV| < 1.5 MeV/c^. The signal and background 
yields (which normalize the pdf 's in the likelihood function) are determined from the m^^- 
Q fit described above. The decay-time distributions for the backgrounds before smearing 
are taken to be: e~*^^DO for random vr background, e~*/^D36 for > 3-body background, 
e~^l'^Dch for jD^ background, and (5(t) for combinatorial background. Parameter r^ig;, 
is obtained from fitting sideband data, and r^^^ is obtained from MC simulation. These 
time distributions are convolved with resolution functions. For random-vr background (and 
the small multi-body / / background), the resolution function used is the same as 
that for signal; for combinatorial background, the resolution function has the same form but 
parameters f['^\^^\ 5'*^'^°™^'', determined from fitting sideband data. 

The fitting procedure is implemented in steps as follows. We first fit the RS signal region 
using a simple background model to obtain a first estimate of signal resolution function 
parameters. We use this resolution function to fit an RS sideband region, which yields 
ft^\^^\ ^-Icomb)^ and r^jg^ for RS background. We then fit the RS signal region with 

these parameters fixed, which yields /x, /ta^jj, 5, S^aii and, as a check, r^o. We subsequently 
use this resolution function to fit the WS sideband region, obtaining /taii™'^\ ^(comb)^ 5'^™™^'', 
and r^3^ for WS background. Finally, we fit the WS signal region, fixing these background 
parameters and those of the signal resolution function; this yields -R^j, x'^, and y' . 

The above fitting procedure has undergone several checks. In MC simulation, background 
parameters obtained from the sideband region fit describe well the background in the signal 
region. The resolution function obtained from data is very similar to that obtained from 
the MC. The lifetime Tno obtained from the RS signal region fit is 415.1 ±1.4 fs, consistent 
with the PDG value |8|; the of the fit projection is 56.0 for 55 bins. We generated 
MC samples having the same size as the data sample, adding the corresponding amount of 
background, and repeated the fitting procedure. For nine sets of (x'^, y') values (spanning the 
ranges [0,2] X 10-2 and [-2,2] x 10"^ respectively), the fit recovers values consistent with 
the generated values. For these sets we also generated ensembles of "toy" MC experiments, 
i.e., without detector simulation, and smeared the time distributions appropriately. Fitting 
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these experiments shows negligible fit bias. 

To this point in the analysis, all optimization of selection criteria was done "blindly," 
i.e., without fitting WS signal events. We now fix the criteria and fit this sample. Four 
fits are done, yielding the results listed in Table HI For the first fit we require that CP be 
conserved. The projection of this fit superimposed on the data is shown in Fig. |21 the 
of the projection is 71.9 for 60 bins. The central value for x'^ is negative (i.e., outside the 
physical region); thus the most-likely value is zero, and we refit the data fixing x'^ = 0. 
The of this fit projection is 73.2 for 60 bins, which is satisfactory. The y' value is ~2(T 
from zero; when we generate MC experiments with this value (and x'^ = 0), we find that the 
probability of obtaining an x'"^ value as negative as what we measure in the data is 8%. For 
the third fit we allow for CPV and fit the ^K^tt~ and K~Tr~^ samples separately. 

For and i?^ we obtain (0.255 i[J:[]i)% and (0.324 ± 0.052)%, respectively. We calculate 
Aj^ and Aj^j and use Eqs. ((21) and (jHl) to solve for x'^ and y'. Finally, for the last fit we 
assume no mixing or CPV and set x'^ = y' = 0; the of this fit projection is 75.6 for 60 
bins, somewhat worse than for the case of mixing. 

TABLE I: Summary of results from the separate likelihood fits. The 95% CL intervals are obtained 
from a frequentist method (see text) and include systematic errors. 



Fit Case 


Parameter 


Fit Result 

(xl0~3) 


95% CL interval 
(xlO-3) 






l-^"^ -1.00 


x'2<0.81 


No CPV 


y' 


or: 4+11.1 
-10.2 


-8.2<2/'<16 






2.87 ±0.37 


2.7<ii£,<4.0 








Rm < 0.42 


No CPV 


y' 


6.0 ±3.3 




x' = (fixed) 




3.43 ± 0.26 








-80 ± 77 


-250<^£,<110 


CPV allowed 




987 ^380 


-991<A^^<1000 








x'2<0.89 




y' 




-30<y'<27 




Pm 




Rj^j < 0.46 


No mixing 
or CPV 


Rjj 


3.81 ± 0.17(stat.)t°;5^(syst.) 



To obtain 95% CL limits on x'^ and y', we use a frequentist method (similar to that 
used in Ref. 01) with Feldman- Cousins ordering P. For points a = (x'^, y'), we generate 
ensembles of toy MC experiments and fit them using the same procedure as that used for 
the data. For each experiment we record the difference in likelihood AL = InL^^g^^ — InL(a), 
where -Z^max is evaluated for x'^ > 0. The locus of points a for which 95% of the ensemble 
has AL less than that of the data is taken as the 95% CL contour. This contour is shown 
in Fig. El 

To allow for CPV, we obtain separate 1 — ^0.05 = 77.6% CL contours for (x'"*"^, y'~^) and 
(x'~^, y'~)- We combine points on the (x'^^, y''^) contour with points on the (x'~^, y'~) 
contour and use these values to solve Eqs. (j21) and (jS)) for x'^ and y'. Because the relative 
sign of x'~^ and x'~ is unknown, there are two solutions (one for each sign); we plot both 
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FIG. 2: The decay-time distribution for WS events satisfying — m^o\ < 22 MeV/c^ and 

I Q — 5.9 MeV| < 1.5 MeV. Superimposed on the data (points with error bars) are projections of 
the decay-time fit. 

in the (x'^, y') plane and take the outermost envelope of points to be the 95% CL contour 
allowing for CPV. This contour has a complicated shape due to the two solutions. Because 
the contour includes the point x'^ = y' = 0, we cannot constrain (p at this CL. 

We evaluate systematic errors by varying parameters used to select and fit the data within 
their uncertainties and recording the new fit values obtained. These values are shifted 
with respect to the original central value Sq. We find the significance m of a shift via the 
formula = — 2[lnL(aj^g^) — lnL(aQ)]/2.30, where the factor 2.30 corresponds to 68.3% 
confidence in two dimensions. We add in quadrature the significances of individual shifts 
to obtain an overall scaling factor ^1 + Y.f^i- We increase the 95% CL statistical contour 
by this factor to include systematic errors. As a check, we generate an ensemble of toy 
MC experiments with aj^c~ 0.006) and fit them to confirm that 68.3% of the ensemble 
satisfies — 2[lnL(^)— InL(ap^Q)] <2.3, where ^ is the central value for an experiment. 

The parameters varied include kaon and pion identification criteria, the allowed of 
the vertex fit, background shape and normalization parameters, and resolution function 
parameters for both signal and combinatorial background. The largest shift in (x'^, y') occurs 
for the D*~^ momentum cut; when this is varied over a significant range, {Ax''^/x''^)^^^^ = 

12%, (Ay7y')m^^ = 10%, and A(-2 InL) = 0.092. The overall scaling factor is ^Jl + E = 

1.08. For the general case allowing for CPV, we scale the D° and contours separately 
before combining. The rescaled 95% CL contours are shown in Fig. |3| the dashed contour 
corresponds to the CP-conserving case and the solid contour to the general case. Projecting 
these contours onto the coordinate axes gives the 95% CL intervals for x'^ and y' listed in 
Table HI _ 

In summary, we have searched for D^-D^ mixing and CP violation using WS ^K^tt~ 
decays. In 90 fb~^ of data we find no evidence for these processes and constrain the mixing 
parameters x'^ and y' and the CP asymmetry parameters Aj^ and Aj^. The limits for x'^ 
and y' are more stringent than previous results; the value for Rj^ (the ratio of DCS to CF 
decays) has smaller uncertainty. 
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FIG. 3: 95% CL regions for (j;'^, y'). The dotted (dashed) contour is statistical (statistical and 
systematic) and corresponds to CP conservation. The dash-dotted (solid) contour is statistical 
(statistical and systematic) and allows for CPV . The open circle represents the most-likely value 
when CP is conserved and x'^ is constrained to be > 0. 
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